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ABSTRACT 

Large galaxy redshift surveys have long been used to constrain cosmological models and structure 
formation scenarios. In particular, the largest structures discovered observationally are thought to 
carry critical information on the amplitude of large-scale density fluctuations or homogeneity of the 
universe, and have often challenged the standard cosmological framework. The Sloan Great Wall 
(SGW) recently found in the Sloan Digital Sky Survey (SDSS) region casts doubt on the concordance 
cosmological model with a cosmological constant (i.e. the flat ACDM model). Here we show that 
the existence of the SGW is perfectly consistent with the ACDM model, a result that only our very 
large cosmological Y-body simulation (the Horizon Run 2, HR2) could supply. In addition, we report 
on the discovery of a void complex in the SDSS much larger than the SGW, and show that such 
size of the largest void is also predicted in the ACDM paradigm. Our results demonstrate that an 
initially homogeneous isotropic universe with primordial Gaussian random phase density fluctuations 
growing in accordance with the General Relativity, can explain the richness and size of the observed 
large-scale structures in the SDSS. Using the HR2 simulation we predict that a future galaxy redshift 
survey about four times deeper or with 3 magnitude fainter limit than the SDSS should reveal a largest 
structure of bright galaxies about twice as big as the SGW. 

Subject headings: cosmology: large-scale structure of the universe - cosmology: observations - galaxies: 
statistics - methods: numerical - methods: observational 



1. INTRODUCTION 

The SGW (Gott et al. 2005) found in the SDSS (Aihara 
Q\ \ et al. 2011) is a thick filamentary structure of galaxies lo- 
l/") ■ cated at a distance of about 300 Mpc from the Earth. Its 
^sO " densest part spans about 200 Mpc, and the whole fila- 
l/"") \ ment projected on a slice appears to be contiguous over 
. the scale of more than 400 Mpc. It is likely to be longer 
^ than the survey size since the structure is cut by the 
\ survey boundaries. The SGW is reminiscent of the CfA 
. Great Wall (Geller & Huchra 1989), that triggered an 
r* 1 intense dispute against the 'Standard' Cold Dark Matter 
(SCDM) model (White et al. 1987) or even the whole 
k> . class of the gravitational instability models of structure 
formation (Park 1990 and references therein). The cos- 
' mological principle of homogeneity and isotropy of the 
universe was also doubted since structures as big as the 
survey size were always found as the survey size was in- 
creased. 

The skepticism was relieved when Park (1990) demon- 
strated that large-scale structures (LSS's) with size up 
to 200 Mpc can appear in the standard CDM cosmolog- 
ical model in surveys like the CfA survey. About 20 
years later the SGW revived the doubt on the cosmolog- 
ical principle and the current standard model (ACDM 
model) since it is even larger than the CfA Wall and 
comparable with the survey size again (Sheth & Diaferio 
2011). 

In this paper we ask whether or not the existence of 
the LSS's observed in the SDSS is consistent with the 
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ACDM model that adopts a homogeneous isotropic uni- 
verse with primordial Gaussian random phase density 
fluctuations as predicted by inflation and a reasonably 
successful galaxy assignment scheme. The comparison of 
real superclusters with simulated ones has been made in 
several previous studies (Einasto et al. 2006; Einasto et 
al. 2007b; Einasto et al. 2007c; Araya-Melo et al. 2009). 
We also make predictions on the properties of the LSS's 
to be observed in the future deeper surveys. Hubble con- 
stant of 72 km s -1 Mpc -1 is used in this paper. 



2. THE SDSS SAMPLE 

To identify the LSS's in the observational sample we 
use the SDSS Main galaxy sample (Strauss et al. 2002), 
which is currently the largest three-dimensional sample 
of galaxies with a high sampling density. A volume- 
limited subsample of 116,877 galaxies with absolute r- 
band magnitude brighter than —21.6 is generated from 
the KIAS value-added catalog (Choi et al. 2010a) , which 
supplements the bright galaxies missing in the SDSS 
sample. The magnitude limit is about 0.6 magnitude 
brighter than the critical magnitude M* of the SDSS 
Main galaxy sample (Choi et al. 2007), and corresponds 
to the sample depth of 689 Mpc. The mean separation 
between galaxies is 12.5 Mpc. The sample is large enough 
to reduce the cosmic variance for the number of SGW- 
like structures, and yet the galaxy number density is high 
enough to trace major LSS's. We calculate the comoving 
distances R of galaxies using the WMAP 5 year cosmo- 
logical parameters (Komatsu et al. 2009) , and the Carte- 
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sian coordinates are calculated as in Park et al. (2007), 

x = — iteinA, y = i?cosAcos?7, z = RcosXsmr), (1) 

where A and rj are the SDSS survey coordinates. 

3. IDENTIFICATION OF THE LARGE-SCALE 
STRUCTURES 

Superclusters of galaxies have been identified by many 
previous studies, which typically use the smoothed lumi- 
nosity density of galaxies and apply a threshold level to 
define structures (Basilakos, Plionis, & Rowan-Robinson 
2001; Einasto et al. 2007a; Luparello et al. 2011; Liiva- 
magi et al. 2012). Here we adopt the Friend-of- Friend 
(FoF) algorithm to identify the high-density LSS's by 
connecting close galaxies because the results convey the 
visual impression well and the method uses effectively 
only one free parameter. 

Before we search for structures in a sample of galaxies, 
we apply the FoF algorithm with the linking length of 
3000km _1 s to the sample to find massive groups. The 
dispersion of each group along the line of sight is forced 
to be equal to that across the line of sight if the former 
is larger than the latter. We then search for LSS's by 
connecting galaxies separated by less than the connec- 
tion length d c . A very small d c results in no LSS's and a 
very large d c gives just one LSS connecting all galaxies. 
We choose to use the critical linking length that results 
in the maximum number of structures (Basilakos 2003). 
In the case of our volume-limited sample we find the crit- 
ical linking length of d c — 7.78 Mpc when the minimum 
number of member galaxies is set to 20. 

The left column of Fig. [1] shows the four richest high- 
density LSS's found in the SDSS sample using the FoF 
method. The richest one among all 873 structures dis- 
covered within our sample is a dense part of the SGW 
containing 822 bright galaxies (the top left panel) . How- 
ever, the largest one is the second richest structure with 
the maximum extent of 226 Mpc containing 601 galax- 
ies. The part of the SGW identified by our FoF method 
is 207 Mpc long. The richest structures are typically fil- 
amentary, and filaments are often looping. Also occa- 
sionally found are superclusters attaching a few spread- 
ing filaments. The relatively poor structures are usually 
single superclusters or doubles connected by a thin fil- 
ament. Each high-density LSSs contains a few super- 
clusters found by several previous studies (Einasto et al. 
2011a; Einasto et al. 2011b; Liivamagi et al. 2012). 

The low-density LSS's are found in a similar fashion. 
We tessellate the SDSS region with cubic pixels, and then 
mark the void core pixels from which there is no or only 
one galaxy within a radius of rbaii- Note that the high- 
density structures are regions where a ball with a radius 
of d c contains two or more galaxies. Therefore, our def- 
inition for void core pixels is consistent with and com- 
plementary to the definition of the high-density struc- 
tures. To identify voids these void pixels contacting their 
faces are grouped together by using the FoF algorithm. 
We then expand each void out to the distance of rbaii 
- d c to account for the buffer region with high-density 
structures, where d c is the linking length used to iden- 
tify the high-density structures. We adopt the critical 
rball = 18.14 Mpc that results in the maximum number 
of independent voids with volume of the core pixel region 
larger than (15 Mpc) 3 . 



On the right of Fig. [T] the three largest-volume void 
complexes found in our SDSS sample are shown. The 
largest one among the 385 voids found is located just 
behind and above (in y and z directions) the SGW (the 
top right panel, see also Einasto et al. 2011c, 2012 for a 
large void at similar location) . Its total volume is as big 
as (157 Mpc) 3 and its maximum extent is 464 Mpc, more 
than twice larger than the dense part of the SGW. The 
largest few void complexes have very complicated topol- 
ogy showing many voids multiply connected by tunnels. 
The relatively smaller low-density structures are usually 
single voids. 

4. TV-BODY COSMOLOGICAL SIMULATION 

Our main concern is to know how often these largest 
structures are expected to be observed in the currently 
popular model of the universe. Taking one more step 
further we ask if the whole distribution of richness or size 
of the observed LSS's is consistent with the prediction of 
the standard model. 

We use a very large cosmological N-body simulation, 
Horizon Run 2 (HR2) (Kim et al. 2011), for the com- 
parison. The simulation evolved 6000 3 particles in a box 
with a side length of 10 Gpc to calculate the gravitational 
evolution of primordial density fluctuations generated in 
accordance with a ACDM model (Komatsu et al. 2009). 
The matter, baryon, and cosmological constant density 
parameters are set to 0.26, 0.044, and 0.74, respectively 
(see Kim et al. 2011). The minimum mass of dark matter 
subhalos identified in the simulation is 5.2 x 10 12 M Q , and 
the mean subhalo separation is 12.5 Mpc, equal to that 
of our SDSS galaxy sample. We assume that each dark 
matter subhalo above the minimum mass contains one 
galaxy. This subhalo-galaxy one-to-one correspondence 
(sHGC) model has been proved to work well in terms of 
one-point function and its local density dependence (Kim 
et al. 2008), two-point function (Kim et al. 2009), and 
also topology (Gott et al. 2009; Choi et al. 2010b). Using 
this galaxy assignment scheme and the idea of abundance 
matching we assume that the subhalos with mass above 
5.2 x 1O 12 M compare with our SDSS galaxies brighter 
than M r = -21.6. 

The total volume of our SDSS sample is (615 Mpc) 3 ef- 
fectively. Our Horizon Run 2 simulation, about 16 times 
larger in linear size, is for the first time large enough 
to capture the large scale power actually present in the 
standard model of cosmology and at the same time has a 
mass resolution high enough to simulate the SDSS Main 
Galaxy Sample. This uniqueness of the simulation en- 
ables us to estimate the statistical likelihood of the LSS's 
found in the observations. 

5. COMPARISON BETWEEN THE OBSERVATIONS 
AND THE ACDM MODEL 

We make 200 SDSS-like surveys of the 'galaxies' in the 
simulation at the present epoch, and analyze the mock 
survey samples in exactly the same way the observational 
data is analyzed. 

For each mock sample the fingers-of-god are identified 
and compressed, and the critical linking length giving the 
maximum number of LSS's is found. We find that the 
mean and standard deviation of the linking lengths from 
the 200 mock SDSS samples are d c = 7.71 ± 0.18 Mpc, 
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Fig. 1. — (Left) The four richest high-density large-scale structures found by the Friend-of- Friend method with the linking length of 
7.78 Mpc. Galaxies belonging to each structure are projected on to the x-y plane of the SDSS survey coordinates (Choi et al. 2010b). L is 
the maximum extent of each structure. (Right) Three largest volume low-density large-scale structures (void complexes). The total volume 
V is calculated by expanding 10.4 Mpc to all directions from the core region shown in the plot to take into account the boundary regions 
with the high-density structures. 



quite close to that of the observational sample. Fig. [2] 
shows the four largest typical high-density LSS's and 
three largest typical low-density structures selected from 
the 200 mock samples. For example, the largest one in 
the hgure is the structure having the approximately me- 
dian maximum extent and also the median richness (or 
volume in the case of low-density structures) among the 
200 largest structures found in the 200 mock samples. 
And the second largest one is the median structure se- 



lected from the 200 second largest ones. 

Out of the 200 mocks 137 samples contain a high- 
density structure richer than the SGW. We also find 
that the largest high-density structure is longer than the 
SGW identified in the same way in 155 cases. Therefore, 
we conclude that structures like the SGW can be easily 
found in surveys like the SDSS in the ACDM universe 
even though the LSSs grew from primordial Gaussian 
fluctuations in a homogeneous isotropic background. On 
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Fig. 2. — The four largest typical high-density LSS's (left) and three largest typical low-density structures (right) selected from the 200 
mock samples. The high-density LSS's show superclusters and filaments connected quite similarly to the observed ones. Likewise, the 
low-density structures show topology of voids very much like the observed void complexes. 



the other hand, none of the mock samples has the 6- 
th richest or largest structure richer or larger than the 
SGW. This means that even though the SGW-like struc- 
tures can be found quite often in a SDSS-like survey, such 
large structures are always one of the top six richest and 
largest structures. The SGW is indeed a rare object, and 
was able to be found because of the large survey volume 
of the SDSS. 

Our conclusion is opposite to that of Sheth & Diaferio 
(2011), who used the extreme value statistics to estimate 



the likelihood of finding a SGW-like object in the SDSS. 
They claimed that the existence of the SGW is a 4cr event 
in the flat ACDM universe with the r.m.s. amplitude of 
density fluctuation in a 11.1 Mpc radius sphere of ag — 
0.8, and is difficult to reconcile with the model. 

To further inspect the consistency between the obser- 
vations and the ACDM model we calculate the distribu- 
tion functions of the richness and size of the LSS's. Fig. [3] 
shows the distributions of the number of structures with 
member galaxies more than N g (open circles in the top 
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Fig. 3. — Richness (upper) and size distributions (bottom) of 
high-density large-scale structures in the observations and simula- 
tions. Cumulative distribution functions of the richness and size of 
the high-density LSS's identified in the SDSS sample (circles). The 
richness N g is the number of member galaxies belonging to each 
LSS. The triangles are the distribution functions when the linking 
length is increased by 10%. The corresponding solid lines and the 
error bars attached are the mean and the standard deviations ob- 
tained from 200 mock surveys in the Horizon Run 2 simulation. 
The red arrows on the left are the locations of the SGW when the 
linking length d c is set to the critical value while the right ones 
are when d c is 10% larger. Note the SGW is the 5th richest and 
largest structure for this larger linking length. 



panel) and with a maximum extent larger than L (open 
circles in the bottom panel). The y-axis $ is the num- 
ber of structures per unit SDSS volume. The mean and 
standard deviation of the cumulative histograms from 
200 mock surveys are shown as lower solid lines and er- 
ror bars. It can be seen that the observed richness and 
size distributions agree astonishingly well with the sim- 
ulation. 

Our ACDM simulation tells us that on average the 
richest high-density LSS in the flat ACDM universe 
in a SDSS volume is expected to contain 957 galaxies 
brighter than M r = —21.6 when the linking length is 
7.78 Mpc, and that the typical size of the largest struc- 
ture is 255 Mpc. These values compare with 822 galaxies 
and 226 Mpc for the SDSS sample. Therefore, the largest 
structures in the observations are actually a little smaller 
than the ACDM expectations both in the richness and 
size . The triangles and upper solid lines are the cumu- 
lative histograms for the observation and the simulation 
when the linking length is increased by 10% above the 
critical value. The plots demonstrate that the agreement 
does not depend on the choice of the linking length. 

Springel et al. (2006) claimed that a SGW-like object 
was found in their Millennium Simulation to support 
their view that LSS does not provide the strongest chal- 
lenge to ACDM. However, considering the fact that the 
matter fluctuation power spectrum of their simulation 
deviated largely from the ACDM theory near the simu- 
lation box scales due to incorrect normalization and sta- 
tistical fluctuations (Springel et al. 2005) and also hav- 
ing only one simulation whose volume is roughly equal 




Fig. 4. — Volume (upper) and size distributions (bottom) of voids 
in the observations and simulations. Cumulative distribution func- 
tions of the volume and size of the void complexes identified from 
the SDSS sample (circles). Two hundred mock surveys in the Hori- 
zon Run 2 simulation are used to calculate the mean (solid lines) 
and standard deviations (error bars). 



to SDSS, it would be difficult to draw such a conclusion 
on the prevalence of a SGW-like structure in the SDSS- 
like surveys in the ACDM universe. Our HR2 is about 
4, 300 times larger than the volume of the SDSS sample 
with the fundamental mode more than 15 larger than the 
depth of SDSS, and accurate statistical comparisons as 
presented here are possible. 

The impressive agreement between the observations 
and the ACDM model is also found for the volume and 
size of void complexes. Fig. |4] shows the number of the 
void complexes in the SDSS sample with volume larger 
than V (open circles in the upper panel) and that with 
the maximum extent larger than L (open circles in the 
bottom panel). The cumulative distribution functions 
agree very well with the mean of the 200 mock surveys 
in the HR2 (solid lines with error bars). 

To find how the size of the largest LSS scales as the 
survey size increases we made 27 non-overlapping mock 
surveys in the HR2 simulation having the SDSS angular 
mask but with the outer boundary located at 2, 767 Mpc 
or redshift 0.8. These surveys are about 4 times deeper or 
3 magnitude fainter than the SDSS. We use these mock 
surveys to correctly account for the survey boundary ef- 
fects on the scaling of the LSS. We find that the largest 
LSS typically has the mean number of galaxies of 2, 480 
and the maximum extent of 430 Mpc. Therefore, the 
maximum richness and size will increase by a factor of 
only about 3 and 1.7, respectively, and the universe will 
look more homogeneous over the scale of the survey size. 

6. CONCLUSIONS 

We identify high-density and low-density LSS's in the 
SDSS to test whether or not the current standard ACDM 
model of the universe can explain the observed struc- 
tures of bright galaxies (Choi et al. 2010b). The LSS's 
used in this comparison are those when the character- 
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istic connection lengths result in the maximum number 
of structures. It is found that the richest high-density 
structure is the dense part of the SGW, which is also the 
second largest structure. The low-density LSS's are typ- 
ically much larger than the high-density counter parts, 
and the largest one is found to be 464 Mpc long. 

The Horizon Run 2 simulation of the ACDM universe 
is used to make a set of mock SDSS surveys. Galaxies 
are assigned to dark matter subhalos assuming that each 
halo contains one galaxy and adopting the abundance 
matching (Kim et al. 2008). LSS's are identified in the 
exactly same way that the SDSS data is analyzed. We 
find that the structures with richness and size similar to 
the SGW are usually one of the richest and/or the largest 
structures in the mock samples. To estimate the statis- 
tical significance of the largest observed LSS's and to 
check the consistency of the properties of the observed 
LSS's with those of the structures found in ACDM we 
measure the distribution functions of richness and max- 
imum extent of LSS's identified in the SDSS and the 
mock samples. It is found that the observed distribution 
functions agree with those of the simulation astonish- 
ingly well. We conclude that both observed high-density 
and low-density LSS's have the richness/volume and size 
distributions consistent with the ACDM universe. This 
agreement between the observations and the theoretical 
predictions should be considered as one of the great suc- 
cesses of the ACDM cosmological model. 

Einasto et al. (2006, 2007b, 2007c) compared prop- 
erties of 2dFGRS and SDSS superclusters with those of 
superclusters identified from the Millennium Run mock 
galaxy catalogue. They showed that the geometric prop- 
erties of real superclusters such as the size, the degree of 
asymmetry and compactness, and the mass of the rich- 
est superclusters are similar to those of simulated ones 
(see also Araya-Melo et al. 2009). It was found that 
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the fraction of such extremely massive and richest su- 
perclusters is too small in the simulated samples when 
compared to the observed samples, and the morphology 
of the richest supercluster in the SGW is not recovered in 
simulation (see also Einasto et al. 2011b). Richness and 
luminosity function of LSSs depend sensitively on three 
things: first, the initial conditions such as the primordial 
power spectrum, second, the galaxy (and their luminos- 
ity) assignment scheme, and third, the LSS identification 
method. Therefore, for a fair comparison between obser- 
vations and simulations, it is very important to use the 
same mass objects and to identify the LSS's using the 
same criteria. Whether or not the largest LSS's of the 
ACDM universe have properties different from the ob- 
served ones remains to be studied further. 

We note in this study that the properties of LSS's de- 
pend sensitively on the initial power spectrum and the 
growth of structures and can be a powerful tool to dis- 
criminate among cosmological models and galaxy forma- 
tion scenarios. We plan to further explore the usefulness 
of using the LSS properties in cosmology in the future 
studies. 
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